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Abstract 

A substantial enhancement is found of the b — ► ccs decay rate due to the QCD interaction 
within the cs pair, which enhancement may amount to 30%, or more. Some general 
features of calculation of the QCD radiative corrections in two first orders are discussed. 



The present theoretical understanding of B decays conspicuously runs into a problem 
with explaining the experimental observation of rather low semileptonic branching ratiotQ , 
which requires^ an enhancement of the non-leptonic decay rate by as much as 20% - 30% 
over the existing theoretical predictions. Therefore a more thorough theoretical study of 
the B decays is imperative. In this note the decay b — > ccs is discussed, whose rate is 
suppressed by the presence of two massive charmed quarks as compared to the dominant 
non-leptonic decay b — > cud. An enhancement of the decay b — *> ccs would somewhat 
relax the problem of low semileptonic branching ratio. However, such enhancement would 
also worsen another possible problem, which is perhaps hinted at by the experiment @, 
that is the problem of a low average charm yield per B decayQ. Nevertheless, within 
the present uncertainty of the measured average charm yield in B decays the data can 
still accommodate a substantial enhancement of the decay b — > ccs. In any case, it is 
important to have the potentially essential effects calculated. Here it is found that the first 
order QCD correction to the total rate of the decay b — > ccs due to the interaction within 
the quark pair cs has relative magnitude 5-c S a s /ir with an unusually large coefficient: 
6cs = 4.46 for m c /mb = 0.3, unlike the case of the decay b — > cud, where the similar 
factor due to the ud pair is equal to one. 

The QCD radiative effects in the decay rate of the b quark are usually analyzed within 
the leading log approximation in ln(m^/m&) or in the next-to-leading log approximation 
(see e.g. |0, ||). However, since, \n(m w /m b ) is not really a sufficiently large parameter, 
the non-logarithmic terms may be quite essential, and it might be more reasonable to 
rely instead on complete calculation of the QCD radiative corrections in the first and 
second orders in a s . The error in the logarithmic terms, induced by such truncation of 
the series, is then not more than about 5% (in the total rate), which is not larger than 
other uncertainties in the calculation, in particular not larger than the non- logarithmic 
terms. 

Both the perturbative and non-perturbative effects in the inclusive decay rates of B 
hadrons are conveniently calculated by representing'^ the decay rate through the absorp- 
tive part of the B self-energy, arising in the second order in the weak interaction. For the 
first-order QCD corrections to the rate of the decay b — > cq 2 qi one thus should consider 
the absorptive part arising from all possible unitary cuts of the graphs of three types 
shown in fig.l. (These graphs ignore the 'penguin' contribution, arising for qi = q 2 , which 

1 Clearly, it is impossible to solve simultaneously both these problems by an enhancement of the decay 
b — > ccs H 
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is to be discussed separately.) In fact, however, it is clear that the graphs of the type 
in fig.lc, i.e. with gluon exchange between the be line and the q 2 qi loop, are vanishing 
because of the color trace over the loop (a gluon can not interact through the loop with 
two colorless W bosons)^. The graphs of the type in fig. la do not involve the interaction 
of the quarks q± and q 2 in the loop with gluons. Therefore the correction arising from 
these graphs can be adapted^ from the old calculations^ of QED corrections to the 
muon decay (see also |T0| , |TT[1). As to the effects of the gluon exchange within the colorless 
q 2 qi loop, these were discussed thus far for both quarks being massless, in which case the 
correction reduces to the familiar factor a s /Ti. This limit is justified for the case of the 
decay b — > cud, but, as will be shown, is misleading for the decay b — > ccs, where the 
charmed quark in the loop has mass, which is not small in comparison with m;, . 
Thus in the first order in a s the rate of the decay b —>■ ccs can be written as 

Tecs ^ ccs ( 1 \$bc ~\~ &cs ^penguin] ) j (1) 

V 7T / 

where is the rate without any QCD corrections, 5b c arises from gluon exchange on 

the be line, 5c S is due to gluon interactions within the cs loop, and S pengu i n is due to effects 
of the penguin type. It is the goal of the present note to calculate the correction factor 

Starting with the relevant term in the weak Lagrangian of the form 

L int = 2 V2 G F V cb V g2gi (c L 7^ b L ) (g x L 7 M g 2 l) (2) 
and parametrizing the spectral density of the current = {^ilI^l) as 

£<°I^(-*><*OZ)|0> = -^-Mq*) Ug, v - q,q v ) + ^-B(q*)q,q u , (3) 

n 07T v ' Q7T 

one can write the total decay rate of b — > cq 2 <?i as 



T c - q2qi = 6r m fe - 8 / {A(q 2 ) q 2 (m 2 b + m 2 c - q 2 ) + 

J(mi+m 2 ) 2 

\ (A(q 2 ) + B{q 2 )) [{m 2 b - m\) 2 - q 2 (m 2 b + m 2 c )}\ \j \{m 2 , q 2 , m 2 c ) dq 2 , (4) 



2 In particular, this explains why there is no first-order correction proportional to a s ln(mi^/mf,): The 
graphs of the types in fig. la (gluon exchange on the be line) and in fig. lb (gluon exchange within the 
loop) contain renormalization of the V — A currents, which is not logarithmic. 
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where X(x, y, z) 
and 



x 2 + y 2 + z 2 



2xy — 2xz — 2yz is the standard kinematical function, 



r 



3G 2 F \V cb V q2qi \ 2 m 5 b 



(5) 



1927T 3 

is the lowest-order parton decay rate with massless quarks in the final state. 

It should be emphasized that the equation (||) is applicable to calculation of only the 
effects associated with the gluon exchanges within the q 2 q\ loop, i.e. the ones, which are 
discussed in this note, and does not include the effects of gluon exchange on the be line 
(like the one shown in fig. la) or the effects of gluon exchange between the loop and the 
be line, which arise starting from order a 2 . In the absence of QCD radiative effects the 
form factors A and B are readily calculable. In the case when q 2 = c and q\ = s, so that 
ni\ — and m 2 = m c , one finds 



An 



2 
3 



1 + 



2q 2 



1 9 \ 2 

q — m c 




(6) 



Using these expressions in eq. 
the 'bare' rate of the decay b - 



|) and integrating over q 2 gives the well known result for 
ccs in eq.(JI|): 



^(bare) 



1 o 1 \ — j — 

,rrib nib 



(7) 



where 



I(x, x) 



Vl - 4:X 2 (1 



Ux 2 - 2x 4 - 12x 6 ) 



+ 24x 4 (1 - a; 4 ) 



In 



1 + Vl-42 2, 
1 - Vl - 



The QCD radiative corrections due to interactions within the cs loop are expressed 
through the radiative corrections to the form factors A and B in the spectral density 
in eq.(||). The calculation of the 0(a s ) corrections to the spectral density with unequal 
masses of quarks has been done in connection with the QCD sum rules both for the 
longitudinal form factor B and for the transversal one v4 0. For the case, 
relevant here, where one of the quarks is massless the result reads as 



A(q 2 ) 
B(q 2 ) 



Ao(q 2 ) 
B (q 2 ) 




(9) 
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with 

fi(z) = — + 2Li(-)+mz In— ^--ln(z-l)+ln— ^ + -ln(z-l) + — -In* (10) 

~r \ Z / Z X ._ Z .L Z Z X 

and 

„ / s 5 1 1 - 1\ /3 1 \ , , . z , , N 

/2<2 > = "2 - 7 - — + (— ) (2 + 2;) '"< 2 - : ) + (731,1 • I") 

where Li(x) = — J X ln(l — t) dt/t is the standard dilogarithm function^. 

The integral in eq.([|) with the radiatively corrected values of the form factors can be 
easily done numerically, thus giving the value of the correction factor Sc S in the equation 
(|I|). The results of such calculation are shown in fig. 2 in terms of the behavior of 5c S as a 
function of m c jm^. In particular, at a 'reference' point m c /mf, = 0.3 one finds 5c S = 4.46, 
which is significantly larger than the analogous correction factor for the b — ► cud: 5ud = 1- 
If one uses the value a s = 0.2, then the discussed correction is about 30%. However, a 
closer inspection of the integral for the correction in eq.([|) shows that the integrand has a 
maximum at q 2 ~ 2m 2 c ~ 0.2ml for the discussed here range of values of m c j nib- Therefore 
the appropriate value of a s can in fact be larger. Naturally, a quantitative clarification of 
this point requires a higher order calculation^. It can be also noted that the enhancement 
of the contribution of relatively low values of q 2 is due to the logarithmic growth of the 
function fi(z) in eq. flToD in the threshold region z — *■ 1, which is a consequence of the 
'hybrid' anomalous dimensionlI3 of the current (cTs). 



We therefore conclude, that the decay b — > ccs is enhanced by about 30% or more 
by the correction proportional to 8-^. To assess the resulting fraction of this decay in the 
total decay rate one should also take into account the corrections with 5 bc and 5 penguin 
and measure the result against, say, the semileptonic mode b — > civ with I = e or I = /i, 
whose rate contains only the QCD correction associated with the be line i.e. with 8®. 
The penguin effect is negative and is about 3% - 5% in mag nitude@. The term 5^ is also 
negative, but due to the charmed quark mass its magnitude is somewhat smaller than 
that of the negative 8®. In effect the term with 5 pengu i n approximately cancels against 

3 The correction to the longitudinal form factor B coincides with the correctionliU for scalar or pseu- 
doscalar density up to an additive constant, corresponding to the normalization of the (pseudo)scalar 
operator in order a s . This constant is 

fixed OJ EJ 

unambiguously for the longitudinal part of the vector 
or axial current. I am thankful to P. Ball for pointing out to me the papers where this point is 

clarified. 

4t 



It has been recently argued 1^ that the natural normalization scale for a s in 5b c is also quite low. 



4 



the difference Sb c — <5j, c m the ratio T C c S /T c i u and the net 0(a s ) correction is dominated by 
the large 5c S - Therefore it is quite possible that a sizeable part of the existing theoretical 
deficiency of nonleptonic decays of B can be eliminated by a 30% or larger enhancement 
of the decay b — > ccs. 

In order to completely quantify the issue of QCD radiative effects in the B decay rates 
and to possibly achieve an accuracy of about 5% in theoretical predictions for the rate of 
each inclusive mode a complete calculation in second order in a s is needed. Though, no 
attempt of such calculation is done in this note, I would like to conclude with a simple 
general remark concerning a calculation of the 0{a 2 s ) corrections to the inclusive decay 
rates of the b quark by the unitary cuts of graphs similar to those in fig. 1, albeit in the 
next order of the QCD perturbation theory Also in that order one can split the graphs 
into few classes. One class is where the gluons are attached only to the quarks on the 
be line. For the dominant decay b — > cud this correction cancels in the ratio T^d/^dv 
Another is where the gluon corrections are fully contained within the q 2 Qi loop, which are 
reduced to the corrections to the form factors A and B in the spectral density (3) . For the 
pair of massless quarks ud this can be read off the corresponding calculation^ for e + e~ 
annihilation into light hadrons. Third class is where one gluon is exchanged on the be line 
and the other within the g 2 <?i 1°°P an d is thus a product of the first-order corrections. 
Finally, because of the color trace over the loop, the gluon exchange between the be line 
and the q 2 qi loop gives a non- vanishing result only when there are two gluons exchanged, 
each starting on the be line and ending on the q 2 q\ loop, an example of such graph is shown 
in fig.3. It is with the latter graphs that the terms proportional to (a s \ia(mw / mb)) 2 and 
to oil hi(mw/mb) are associated and which may also contain large non-logarithmic terms. 

I am thankful to M. Shifman and A. Vainshtein for discussions of the B decays and 
to M. Danilov for a discussion of the experimental situation with the average charm yield 
in B decays. This work is supported, in part, by the DOE grant DE-AC02-83ER40105. 

When this paper was finished, a revised version of Ref.|| appeared, where a similar 
estimate of 30% enhancement of T C c S was found. I thank P. Ball and V.M. Braun for point- 
ing out to me their revised estimate. As discussed in this note, the actual enhancement 
can in fact be larger due to a larger value of a s at the relatively low relevant invariant 
mass of the cs pair. 
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Figure 1: Three types of graphs, whose unitary cuts describe the first QCD radiative 
corrections to the inclusive decay rate b — > cq 2 Qi- The small filled circles represent the W 
boson propagators and the dashed lines correspond to gluons. The gluon vertices can be 
anywhere on the be line (a), quark lines in the loop (b), or one vertex anywhere on the be 
line and the other vertex on either line in the loop (c). 
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Figure 2: The correction factor 5c S in eq.(l), arising from the gluon exchange within the 
cs loop, vs. the mass ratio m c /m{,. The range of m c /mt, shown well covers the ratio of 
the actual quark masses with the existing uncertainty. 
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Figure 3: The only type of graphs in the second order in a s with gluon exchange between 
the be line and the q 2 qi loop, which gives contribution to the inclusive decay rate. Each 
of the two gluons should start anywhere on the be line and end on either of the quark 
lines in the loop. 
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